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STATISTICAL  ANALYSIS  OF  THE  STABILITY  OF  IMPERFECT  CYLINDRICAL 


SHELLS 

B.  P.  Hakarov  (Moscow) 

This  work  examine'-,  the  problem  of  stability  of  a thin  cylin- 
drical shell,  the  middle  surface  of  which  has  random  initial  devia- 
tions from  the  ideal  shape.  Tor  a solution,  we  used  the  method  of 
V.  V.  Bolotin  [1],  For  the  beginning,  we  use  statistical  character- 
istics of  initial  errors  obtained  according  to  experimental  data  in 
work  [2],  lie  arc  trying  to  construct  a solution  to  the  stochastic 
problem  of  stability  of  a cylindrical  shell  with  a random  calculated 
number  of  random  parameters  which  characterize  the  deviations  of  the 
middle  surface  from  the  ideal  shape,  de  selected  for  these  parame- 
ters coefficients  of  a double  Fourier  series,  with  the  aid  of  which 
the  errors  of  initial  errors  arc  presented. 
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Lot  us  aasune  that  deviations  of  the  middle  surface  of  the 
Shell  f«on  the  ideal  cvlindrical  shape  can  be  presented  In  the  form 
o"  a "our ie”  series: 

*M*i.  *.)  - 2*<$>.in-=p-  + 2 -^fLsin^L  + 

m m mi  h n 

+ 2 2*«>«in-2p-cos4p.  (5) 

m n n 

i.ere,  fm «(1,i  W*  - coefficients  of  breakdown  whicu  '..ill 

subsequently  be  examined  as  random. 

'e  will  search  ’’or  a solution  in  the  form  of  a series  similar 
to  ( j ) and  satisfies  the  boundary  conditions  of  the  problei: 

(*i.  *1)  - S/S’-ln  + 2 2/mUin  ^J-sin  42-  + 

m m n L,  H 

+ 2 2/2  Uin^S-co.4?- 

m n ^ 

’or  coefficients  v:e  obtain  expression 


(C) 


(9) 


-s*nr  the  solution  in  (G),  we  can  determine  the  prccrltical 
stress  in  the  shell,  Thus,  for  example,  the  surroundinr  stresses  in 
the  middle  surface  will  be  equal  to 

nr  - - 4[2^#,,in  nr1  + 22 /*»»  + 

+ 22 /a  (srnbF»in  -n?- co>  -nr] 

Let  us  set  up  a problem  on  the  stability  of  solution  (G) . The 

equation  of  neutral  eouillbrium  will  be  obtained  by  linearization 
near  the  undisturbed  equilibrium,  here,  as  this  is  done  normally  in 
the  theory  of  elastic  stability,  we  will  consider  the  system  "rigid". 
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disreparllnr  the  precritical  noveaent:;  and  initial  deviations  in 
comparison  with  the  deviations  from  undisturbed  equilibrium  w* . As 
experimental  data  show,  por  shells  which  do  not  have  specially  cre- 
ated defects,  this  linearization  is  fully  ripdit.  The  equation  for 
deviations  w*  take  the  form: 

2>AAAAir.  f-  -jp-  = AA  (10) 

Precritical  surrounding  stress  in  the  middle  surface,  deter- 
mined according  to  formula  (9),  enter  into  the  r ip.ht  half  of  eciua- 
tion  (10).  "'or  stresses  il,„0  v;e  have  expression: 

Uh 

Nu  = — N+  , Nn  *=  , Nlt  *=  Nn  =»  — ( 11 ) 

Let  us  note  that  a similar  approach  was  earlier  used  by  S.  .1. 
Kahn  [3]  for  an  analysis  of  the  stability  of  shells  which  have  axi- 
symmetrical  errors. 

An  exact  analysis  of  the  equation  of  neutral  enuilibrium  (10) 
is  difficult  due  to  the  presence  in  the  riqht  half  of  multiplier 
. Therefore,  for  determininr  the  critical  forces,  we  will  use  a 
variation  of  the  Bubnov-Qalerkin  method. 

Let  us  find  function  w*  in  the  form  of  an  expression  which  sa- 
tisfies conditions  on  the  edres  of  the  shell: 

•»  f%  sin  pjp-  ain  -—jf-  + /*  sin  -^p  coa  pjjr  (12) 

where  1 and  k - wave  numbers  which  characterize  the  form  of  loss  of 
stability.  Expression  (12)  was  selected  with  consideration  of  the 
shift  of  phase  in  the  surrounding  direction. 
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Lot  U3  p 1 n c c (12)  into  equation  (10)  and  demand  it  so  that  it 
will  Le  satisfied  in  terns  of  the  Bubncv-Galerkin  method.  Setting 
to  zero  in  order  to  determine  the  systems  of  equations  with  respect 
tn  fj  and  f , we  obtain  a relationship  which  connects  the  critical 
value  of  the  load  parameter  with  the  coefficients  of  breakdown  of 
the  ''unction  of  initial  errors 
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Coefficients  an  and  depend  on  the  wave  numbers, 
delationship  (13)  can  be  significantly  simplified  if  we  consi- 
der the  followinr  inequality: 

T»».  ^ 1 (10) 

In  this  case.  Tor  v,  we  obtain  an  explicity  expression 

* - + Xt  + y x? + xt») 

where 

* f — • X>-|-^V  (IT) 

Let  us  note  that  in  the  expression  for  v,  the  wave  numbers  ,) 
and  k remained  random,  ^or  their  determination  in  the  future,  we 
must  minimize  v. 

Let  us  employ  the  aonroximation  relationship  (16)  for  studying 
the  statistical  tie  between  critical  stresses  and  random  deviations 
from  the  ideal  3haoe. 


he  will  treat  breakdown  In  (p)  as  a spectral  presentation  of 
randon  function  v?q.  Then  the  set  of  values  £m<°\  &«n.  lHh  will  describe 
the  randon  spectrum  for  Wq.  Let  us  assume  that  these  coefficients 
forn  the  system  of  normal  random  values.  In  this  case,  the  probabi- 
lity functions  for  coefficients  Xg,  , and  X^  are  compositions  of 
normal  laws.  Tor  their  determination,  it  is  sufficient  to  have  a 
correlation  matrix  of  the  spectrum  of  initial  errors. 

' 'e  must  determine  the  mathematical  expectation  of  critical  load. 
Averaging  relationship  (16),  v:e  obtain  in  the  first  approximation 

<v>  * vj*/(l  + <X0>  + </X1*  + Xl»»  (18) 


The  mathematical  expectations  in  the  right  part  arc  expressed 
through  the  characteristics  of  the  coefficients  from  scries  (8)  in 
the  '’oiiowing  manner: 


<Xt>  = 2 <ug»>  , <*»>  = 2 <SJK> 

ftx  mn«  m mn* 

+ X,»  > eXP  {~Ii  — v>ur  * 


(19) 


Dispersions  cr-^'~  and  and  also  coefficient  of  correlation 

v/ ’ * represent  linear  combinations  from  the  elements  of  the  correla- 
tion matrix  of  the  spectrum  of  initial  errors. 

A particularly  simple  result  is  obtained  in  the  case  where 
function  v.'q  is  uniform  in  the  circular  direction.  The  expression 
for  [v]  takes  the  form: 


<*>- *W(1 + <*•>  + V*  < I X,|» 
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In  the  final  st are,  wo  must  minimize  value  [v]  according  to 
the  wave  numbers  ,1  and  , aceom.plinhim’  a computation  for  the  v.’ave 
ranpe  which  correspond:',  to  the  approximate  shape  of  loss  Ir  stabi- 
lity. 

■Relationships  like  (lb)  can  be  used  also  for  determining  the 
probability  density  of  critical  load,  i’ol  lowing'  v:ork  [1],  we  obtain 
the  following’  expression  for  p(v): 

P (v)  = « 2 P [g  (v.  Xt,Xt)Xu  Xt ] |~— 1 dXidXt  (21) 

where  c - normalizinr  multiplier,  and  function  c(v,  X-,,  X-j ) enuals 

* (V.  Xt,  Xt)  = Vjh/V  - 1 -Yx f=T| 

Let  us  briefly  stop  at  sore  numerical  results.  Lot  us  examine, 
for  example,  the  series  op  cylindrical  shells,  for  which  similar 
experinenatal  investigations  wore  conducted  In  work  [2],  nominal 
values.  op  parameters  were  as  follows: 

h - 0.018  cjk,  R « 10  c*.  L = 30  cm,  E = 2.1  10*  kq/C™* 

The  computations,  conducted  with  various  values  of  wave  number 
provided  the  followinr  results: 

* — 3 4 5 0 7 8 

<v>  — 1.12  0.92  0.85  0.42  0.44  0.31  (Xj/k—  1) 

<v>  — 1.10  1.21  0.95  0.81  0.70  0.58  (Xj/*  =»  1,5) 

The  smallest  value  [v]  corresponds  to  k=o,  X^/k=l,  According- 

to  the  data  in  the  experiment,  [v#]*0,23. 

The  difference  between  the  theoretical  and  experimental  values 
was  30/1,  which  attests  to  the  satisfactory  congruence  of  results. 

'•’or  constructing^  the  distribution  of  critical  forces,  v;e  can 
use  relationship  (21).  Integration  In  the  form  of  (21)  is  easiest 
to  accomplish  numerically. 
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